MONODROMY OF HYPERSURFACE SINGULARITIES 



MATHIAS SCHULZE 

Abstract. We describe algorithmic methods for the Gauss-Manm 
connection of an isolated hypersurface singularity based on the mi- 
crofocal structure of the Brieskorn lattice. They lead to algorithms 
for computing invariants like the monodromy, the spectrum, the 
spectral pairs, and M. Saito's matrices A and A\. These algo- 
rithms use a normal form algorithm for the Brieskorn lattice, stan- 
dard basis methods for power series rings, and univariate factoriza- 
tion. We give a detailed description of the algorithm to compute 
the monodromy. 
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1. Introduction 

We consider a germ of a holomorphic map / : (C n+1 , 0) > (C, 0) 

with isolated critical point and Milnor number /i. J. Milnor [Mil68] 
first studied this situation by differential geometry. The regular fibres 
of a good representative over a punctured disc form a fibre bundle 
with fibres of homotopy type of a bouquet of fi n-spheres. The coho- 
mology of the fibres form a flat vector bundle and there is an associated 
flat connection on the corresponding sheaf of holomorphic sections, the 
Gauss-Manin connection. Moreover, there is a monodromy represen- 
tation of the fundamental group of the base in the cohomology of the 
general fibre. A counterclockwise generator acts via the monodromy 
which is an automorphism defined over the integers. Using the De 
Rham isomorphism, the cohomology of the fibres can be described in 
terms of holomorphic differential forms. 

E. Brieskorn [Bri70] first gave an algorithm to compute the com- 
plex monodromy based on this idea. The original Brieskorn algorithm 

l 




was first implemented by P.F.M. Nacken [Nac90] in the computer al- 
gebra system Maple V. A later implementation [Sch99, SchOlc] by 
the author in the computer algebra system Singular [GPS01] using 
standard basis methods turned out to be more powerful. 

In section 2 and 3, we briefly introduce the monodromy, the Gauss- 
Manin connection, the V-filtration, and the Brieskorn lattice, and sum- 
marize properties and results which are relevant for the computation 
of the monodromy. This will lead us to the notions of saturation and 
resonnance of lattices in section 4. In section 5, we describe algorith- 
mic methods for the Gauss-Manin connection based on the microlocal 
structure of the Brieskorn lattice and the Fourier-Laplace transform 
[SchOO, SS01]. We use standard basis methods, univariate factoriza- 
tion, and a normal form algorithm for the microlocal structure of the 
Brieskorn lattice, the latter of which is not published yet. These meth- 
ods lead to algorithms [Sch02, SchOla, SchOlb] to compute Hodge- 
theoretic invariants like the spectrum, the spectral pairs, and M. Saito's 
matrices A and A 1 [Sai89]. We describe an algorithm to compute the 
complex monodromy based on the above ideas. This is also implemeted 
in Singular [SchOlb] and is much faster then the original Brieskorn 
algorithm [SchOlc]. In section 6, we explain how to use the Singular 
implementation by computing an example with a Jordan block of size 
3 x 3 of the monodromy. 



MONODROMY OF HYPERSURFACE SINGULARITIES 3 
2. MONODROMY AND GAUSS-MANIN CONNECTION 

Let / : (C n+1 , 0) > (C, 0) be an isolated hypersurface singu- 
larity. We choose local coordinates x = (x , . . . ,x n ) at G C n+1 and 
t at G C and set d := (<9 , . . . , d n ) where dj := d Xj . We denote by 

u:=dim c C{x}/(d(/)> <oo 

the Milnor number of /. We choose an (n + l)-ball B centered at 
G C n+1 and a disc T centered at G C such that 

Bn/ _1 (0) =: X— f -^T 

is a Milnor representative [Mil68]. Let T := T\{0}^^ T be the 
inclusion of the punctured disc. Then the restriction 

x\f-\o) =: X' — 

is a "io 00 fibre bundle, the Milnor fibration. The fibres X t := f~ l (T) 
have homotopy type of a bouquet of u n-sheres and hence the coho- 
mology of the general fibre is given by 



a\x t ,z) 



W, k = n, 
0, else. 



The n-th cohomology bundle 

H n := R n f*C x , = U Hn ( X ^ C ) 
teT' 

is a locally constant sheaf. The Gauss-Manin connection is the 

associated flat connection 

on the sheaf of holomorphic sections := @t' ®g t , H n which is 

defined by 

V(g <S> v) := dg <g> v, g G & T ', v G H n . 

We denote by d t :— Vg t ■ ^ n > ,ffl n its covariant derivative with 

respect to d t . Lifting paths in T' along flat sections defines an action 

7Ti (T", i) ► Aut (R n (X t , Z)) 

of the fundamental group 7Ti(T", £) on the n-th cohomology of the gen- 
eral fibre. A counterclockwise generator of iri(T',t) acts by the mon- 
odromy 

M G Aut(H n (X t ,Z)). 

Theorem 2.1 (Monodromy Theorem). TTie eigenvalues of the mon- 
odromy are roots of unity and its Jordan blocks are of size at most 
(n + 1) x (n + 1). 
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Let u : T°° > T , u(t) :— exp(27rir), be the universal covering of 

T' where r is a coordinate on T°°. The canonical Milnor fibre is 
defined to be the pullback 

X°° :=X' x T , T°° 

to the universal covering. The natural maps X u ^ = X^°c y x°° 

are homotopy equivalences. We consider A e H ri (X°°,C) as a global 
flat multivalued section A(t) in Jt? n . Note that 

d t A(t) = 0, 

M(A)(t) = A(t + 1). 

3. V- Filtration and Brieskorn Lattice 

Let M = M S M U be the decompostion of the monodromy into semisim- 
ple and unipotent part and set 

N := \ogM u . 

Note that N n+1 = by the monodromy theorem. Let 

H n (X°°, C) = H n (X°°, <D) A , H n (X°°, C) A := ker(M s - A), 

A 

be the decomposition into generalized eigenspaces of M and 

Ma := M|h"(x°°,c) a - 
For A e R n (X°°, C) A , A = exp(-27ria), a e Q, we set 

MM*) :=t a exp^-^-Aog?! A(t). 

Then ip a (A) is monodromy invariant and hence a global section in J(? 71 . 
Definition 3.1. We call the C{t}[9t]-module 

:= {{i^ a {A)) \a GQ,AG H n (X oo ,C) cxp( _ 2 ™ ) >^ C (i*3T) 
the Gauss-Manin connection of /. 
For all a e Z, the map 

^:H"(X-,C)a c ^ 

is an inclusion with image C a := imip a . The following lemma shows 
the correspondence between the monodromy action on H n (X°°, C) and 
the C{t}[<9 t ]-module structure on <S via the maps ip a . 

Lemma 3.2. 

(1) to^) a = ip a+1 

(2) d t o^ a ^^ a _ 1 o(a-£ i ) 

(3) (^-a)o^ = ^o(-^) 

(4) exp(-27rit9i) o ip a — ip a o Ma . 

(5) C a = kei(td t - a) n+l 
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(6) t : C a > C a+1 is bijective. 

(7) d t : C a > C^ 1 i s bijective for «^0. 

The generalized eigenspaces C a of the operator td t define a filtration 
on . 

Definition 3.3. The V-filtration V on is the decreasing filtration 
by C{t}-modules 

V a := < C{t}C' 3 , V >a := 

a</3 a</3 

Note that V a /V >a = C a . There is not only the C{t}-structure on 
. For a > — 1, the action of d t _1 on V a extends to a structure over 

a power series ring C{{9 t -1 }}. This structure is the key to powerful 

algorithms. 

Definition 3.4. The ring of microdifferential operators with con- 
stant coefficients is defined by 

^ fe>0 ' ' 



k>0 



Note that C{{^ -1 }} is a discrete valuation ring and t a <C{t}, a G Q, 
is a free C{{9 t _1 }}-module of rank 1. Together with lemma 3.2, this 
implies the following proposition. 

Proposition 3.5. 

(1) For all a e Q, V a is a free C{t}-module of rank fi. 

(2) & is a ii- dimensional C{t}[t~ 1 ] -vector space. 

(3) For a > —1, V a is a free C{{(9 t -1 }} -module of rank \i. 

There is a lattice in on which the action of d t can be computed. 
Definition 3.6. 

^"^ng/d/Ado- 1 . 

is called the Brieskorn lattice. 

The Brieskorn lattice is embedded in 5f and is a €{t}- and C{{<9f _1 }}- 
lattice. There is an explicit formula for the action of d t in terms of 
differential forms. 

We summarize these well known properties of the Brieskorn lattice 
in the following theorem. 

Theorem 3.7. [Seb70, Bri70, Pha77, Mal74] 

(1) Jf" is a free <C{t}-module of rank ji. 

(2) s: ^"C >Sf ,*(["])(*) 

(3) d t s{[dfA V ])=s{[d V }) 

(4) Jt?" is a free C{{<9 t -1 }} -module of rank /i. 

(5) V- 1 D M"' D V n ~ l 
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The action of d t may have a pole of order up to n + 1 on Jtf". Since 
the monodromy is related to the action of td t , we consider ^-invariant 
lattices in the next section. 

4. Saturation and Resonance 

Since C{t} is a discrete valuation ring, for any two C{t}-lattices 
if, if' C ^, there is a k G Z such that t fc i^ C if'. Hence, for any 
<D{t}-lattice _Sf, 

\/ a2 c =2? c V >ai 

for some «i, a 2 G Q. 

Definition 4.1. Let if C be a C{t}-lattice. IHdjif C JSf then if is 
called saturated and the induced endomorphism £<9 t G End c (if /tif) 
is called the residue of if. If t<9j has non-zero integer differences of 
eigenvalues then if is called resonant. 

Let if C # be a <D{t}-lattice with \/° 2 C JSf C V >ai for some 
Oii,a 2 G Q. By the Leibnitz rule and since V ai is saturated, 

ifo := if ) iffc+i '■= iffc + tdt<S?k, 
defines an increasing sequence of C{t}-lattices 

rcy c^c--c v ai . 

Since V ai is noetherian, this sequence is stationary and 

i-oo • J i-fc- 

fc>0 

is a saturated C{t}-lattice. 

Definition 4.2. ifoo is called the saturation of if. 

The following proposition is not difficult to prove using lemma 3.2. 
Proposition 4.3. [GL73] if M _i = if^ 

If if is saturated then 

if = ( i^nr) ©\/ Q2 . 

ai<o<a2 

Together with lemma 3.2, this implies the following proposition. 

Proposition 4.4. Let if C fre a saturated <C{t} -lattice with residue 
W t E Ead c (Sf/tSf). 

(1) TTie eigenvalues ofexp(—2nitdt) are the eigenvalues of the mon- 
odromy. 

(2) // if is non-resonant then exp(—2nitdt) is conjugate to the 
monodromy. 
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5. Microlocal Structure and Algorithms 

We abbreviate Q* := Q* X o- The microlocal structure extends to a 
C{{s}}[(9 s ]-module structure by the Fourier-Laplace transform 

s := d t ~\ d s := d t H = s-H. 

Since 

[d s , s] = [d t 2 t, or 1 } = d t hdf l - d t t = 1, 

V a for a > -1 and M"' are C{{s}}[s 2 <9 s ]-modules. Note that t = s 2 d s 
is a C{{s}}-derivation. Since 

d t t = s~H = sd s , 

the saturation Jzfoo of a C{t}- and C{{s}}-lattice «Sf is a saturated C{t}- 
and <D{{s}}-lattice. Note that this holds for ". 

By the finite determinacy theorem, we may assume that / e C[i] is 
a polynomial. A C-basis of 

JfC'/sM"' = tt n+1 /df A fi n = €{x}/{d(f)) 

represents a C{{s}}-basis of Jif". If g is a standard basis of (d(f)) 
with respect to a local monomial ordering then the monomials which 
are not contained in the leading ideal (leadg) = (lead(<9(/))) form a 
monomial C-basis of <C{x}/ (df). Hence, one can compute a monomial 

: J of JT". 
m"/ 

We define the m-matrix A = A(s) = ^fc> ^Us fc °^ * ky 

tm =: Am. 

Then the m-basis representation of t on J^" is given by 

tgm = (gA + s 2 d s (g))m. 

If [/ is a C{{s}}[s _1 ]-basis transformation and A' the [/m-matrix of t 
then 

A' = (UA + s 2 d s {U))U- 1 

is the basis transformation formula for U. 

Let resp. J^" be the (rr)-adic resp. (s)-adic completion of Q resp. 

Jtf". The isomorphism da; : C[x] - > f2 n+1 induces an isomorphism 

and M"' is a differential deformation of the Jacobian algebra C[x]/ (d(f)) 
in the sense of [SchOld]. Using the normal form algorithm in 

[SchOld], one can compute any K-jet A< K = A< K (s) = J2k=o^kS k 
of A. 
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The m-basis representation H'^ of the saturation can be com- 
puted recursively by 

H» := €{{s}r, 
HZ +1 :=H»+{s- 1 H»A< k + 8d.HZ). 

Note that only finite jets of A are involved. We use a local monomial 
degree ordering. By computing a normal form of H% +1 with respect to 
H'l, one can check when the sequence 

Hi C E'l C H'J, C • • • 

becomes stationary and compute generators of the saturation H'^ of 
Hg. 

: of 

H 'L is a C{{s}}-basis. The Mm-matrix A' = A'(s) = T, k >o A kS k of t 
is defined by 

MA + s 2 d s M : A'M. 

We set 

5(M) := max{ord(m}) - ord(mf) \m) ^ ^ mf } < /i - 1 
such that 

(MA<^ +5(M) + s 2 «9 s M)^ =: A'<^M. 

for any K > 0. Note that only finite jets of A are involved. Hence, one 
can compute any K-jet A' <K = A' <K {s) = J2k=o^'k sk °f Note that 
the Mm-basis representation of t on is given by 

tgMm = (gA! + s 2 d s (g))Mm. 

Hence, the Mm-basis representation of d t t = s~H on Jf?^/ s K Jf?£ is 
given by 

d t tgMm = (s^gA^ K + sd s {gj)Mm 

and s~ 1 A' <1 = A[ is the Mm-basis representation of the residue of 
J%£. Note that the eigenvalues of A[ are rational by the monodromy 
theorem and can be computed using univariate factorisation. If A[ is 
non- resonant then exp(— 2-rriA^) is a monodromy matrix. 

Otherwise, we proceed as follows. Let S(A') > be the maximal 
integer difference of A[. First we compute A' 1+s , A ,-> from A 1+ s(m)+s{A') 
as before. After a C-linear coordinate transformation, we may assume 
that 



A' 



A' 1 ' 1 A' 1 ' 2 
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with A'\' z = 0, A'f 1 = 0, A' Q = 0, the eigenvalues of A' ' are minimal 
in their class modulo Z, and the eigenvalues of A' 2 ' 2 are non-minimal in 
their class modulo Z. Then the C{{s}}[s _1 ]-coordinate transformation 



U 



s 
1 



gives 



Note that = 0, 5{A") < 5(A') - 1, and that depends only on 
A' <K+1 . With <5(A') of these transformations we decrease S(A') to zero 
such that exp(— 27rL45.) is a monodromy matrix as before. 

The above methods lead to algorithms [Sch02, SchOla, SchOlb] to 
compute Hodge-theoretic invariants like the spectrum, the spectral 
pairs, and M. Saito's matrices A and Ai [Sai89]. 



6. Examples 

The algorithms described above is implemented in the computer alge- 
bra system SINGULAR [GPS01] in the library gaussman. lib [SchOla]. 
We use this implementation to compute an example. 

First, we have to load the library: 

> LIB "gaussman. lib" ; 

Then we define the ring R := Q[x, y, 2]( x , a ,z) and the polynomial / = 
x 2 y 2 z 2 + x 7 + y 7 + z 7 e R: 

> ring R=0, (x,y,z) ,ds; 

> poly f=x2y2z2+x7+y7+z7; 

Finally, we compute the Jordan data of the monodromy of the singu- 
larity defined by / at the origin. 

> spprint (monodromy (f) ) ; 

((1/2,1) ,18) , ((1/2,3) , 1) , ((9/14,1) ,15) , ((9/14,2) ,3) , 
((11/14,1), 15), ((11/14, 2), 3), ((6/7,1), 3), ((13/14,1), 15), 
((13/14,2) ,3) , ((1,2) ,1) , ((15/14,1) ,15) , ((15/14,2) ,3) , 
((8/7,1) ,3) , ((17/14,1) ,15) , ((17/14,2) ,3) , ((9/7,1) ,3) , 
((19/14,1) ,15) , ((19/14,2) ,3) , ((10/7, 1) ,3) , ((11/7, 1) ,3) , 
((12/7,1), 3) 

The computation takes about 2 minutes in a Pentium III 800. A 
Jordan block of the monodromy of size s with eigenvalue exp(— 27rio;) 
occuring with multiplicity m is denoted by ((a, s),m). Note that there 
is a Jordan block of size 3 with eigenvalue —1 which is the maximum 
possible size. 
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